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We show that charged-particles decaying in the early Universe can induce a scale-dependent or 
'running' spectral index in the small-scale linear and nonlinear matter power spectrum and discuss 
examples of this effect in minimal supersymmetric models in which the lightest neutralino is a 
viable cold-dark-matter candidate. We find configurations in which the neutralino relic density is 
set by coannihilations with a long-lived stau, and the late decay of staus partially suppresses the 
linear matter power spectrum. Nonlinear evolution on small scales then causes the modified linear 
power spectrum to evolve to a nonlinear power spectrum similar (but different in detail) to models 
parametrized by a constant running a s — dn s /d\nk by redshifts of 2 to 4. Thus, Lyman-a forest 
observations, which probe the matter power spectrum at these redshifts, might not discriminate 
between the two effects. However, a measurement of the angular power spectrum of primordial 
2f-cm radiation from redshift z ~ 30-200 might distinguish between this charged-decay model and 
a primordial running spectral index. The direct production of a long-lived charged particle at future 
colliders is a dramatic prediction of this model. 

PACS numbers: 98.80.Cq, 98.80.Es, 95.35+d, 12.60.Jv, 14.80.Ly 



I. INTRODUCTION 

While recent cosmological observations provide con- 
vincing evidence that nonbaryonic dark matter exists 
we do not know the detailed particle properties of 
the dark matter, nor the particle spectrum of the dark 
sector. There has been considerable phenomenological 
effort towards placing model-independent limits on the 
possible interactions of the lightest dark-matter parti- 
cle (LDP) 1 in an attempt to try and identify candidates 
within detailed particle theories or rule out particular 
candidate theories. For instance, models with stable 
charged dark matter have been ruled out Q , while signif- 
icant constraints have been made to dark-matter models 
with strong interactions , self- interactions 0, [j| , and 
a millicharge 0, Q ■ Recently, it was shown that a neu- 
tral dark-matter particle with a relatively large electric or 
magnetic dipole moment remains a phenomenologically 
viable candidate 

Concurrent with this phenomenological effort, theo- 
rists have taken to considering physics beyond the stan- 
dard model in search of a consistent framework for viable 
dark-matter candidates. The leading candidates, those 
that produce the correct relic abundance and appear in 
minimal extensions of the standard model (often for inde- 
pendent reasons), are the axion 9] and weakly interact- 
ing massive particles (WIMPs), such as the neutralino, 
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the lightest mass eigenstate from the superposition of the 
supersymmetric partners of the C/(l) and SU(2) neutral 
gauge bosons and of the neutral Higgs bosons PJJ, El • 
However, other viable candidates have also been consid- 
ered recently, such as gravitinos or Kaluza-Klein gravi- 
tons produced through the late decay of WIMPs [l2.[l3| . 
These latter candidates are an interesting possibility be- 
cause the constraints to the interactions of the LDP do 
not apply to the next-to-lightest dark-matter particle 
(NLDP), and the decay of the NLDP to the LDP at 
early times may produce interesting cosmological effects; 
for instance, the reprocessing of the light-element abun- 
dances formed during big-bang nucleosynthesis [bH Il5| , 
or if the NLDP is charged, the suppression of the matter 
power spectrum on small scales and thus a reduction in 
the expected number of dwarf galaxies 0] . 

In this paper we describe another effect charged 
NLDPs could have on the matter power spectrum. If 
all of the present-day dark matter is produced through 
the late decay of charged NLDPs, then, as discussed in 
Ref. 0|, the effect is to essentially cut off the matter 
power spectrum on scales that enter the horizon before 
the NLDP decays. However, if only a fraction of the 
present-day dark matter is produced through the late 
decay of charged NLDPs, the matter power spectrum is 
suppressed on small scales only by a factor (I — /^) 2 . This 
induces a scale-dependent spectral index for wavenum- 
bers that enter the horizon when the age of the Universe 
is equal to the lifetime of the charged particles. What 
we show below is that, for certain combinations of 
and of the lifetime of the charged particle r, this sup- 
pression modifies the nonlinear power spectrum in a way 
similar (but different in detail) to the effect of a con- 
stant a s = dn s /d\nk ^ 0. Although these effects are 
different, constraints based on observations that probe 
the nonlinear power spectrum at redshifts of 2 to 4, such 



2 



as measurements of the Lyman-a forest, might confuse 
a running index with the effect we describe here even if 
parametrized in terms of a constant a s . This has signif- 
icant implications for the interpretation of the detection 
of a large running of the spectral index as a constraint 
on simple single-field inflationary models. The detection 
of a unexpectedly large spectral running in future obser- 
vations could instead be revealing properties of the dark- 
matter particle spectrum in conjunction with a more con- 
ventional model of inflation. We note that the Sloan 
Digital Sky Survey Lyman-a data |l7j has significantly 
improved the limits to constant-a s models compared to 
previous measurements alone [H, Ha Ell- A detailed study 
of the (t,/^) parameter space using these and other cos- 
mological data would also provide interesting limits to 
the models we discuss here. 

While, even with future Lyman-a data, it may be dif- 
ficult to discriminate the effect of a constant running of 
the spectral index from a scale-dependent spectral index 
due to a charged NLDP, other observations may never- 
theless discriminate between the two scenarios. Future 
measurements of the power spectrum of neutral hydro- 
gen through the 21cm-line might probe the linear matter 
power spectrum in exquisite detail over the redshift range 
z « 30 — 200 at comoving scales less than 1 Mpc and per- 
haps as small as 0.01 Mpc [23|; such a measurement could 
distinguish between the charged-particle decay scenario 
we describe here and other modifications to the primori- 
dal power spectrum. If, as in some models we discuss 
below, the mass of these particles is in reach of future 
particle colliders the signature of this scenario would be 
spectacular and unmistakable — the production of very 
long-lived charged particles that slowly decay to stable 
dark matter. 

Although we describe the cosmological side of our cal- 
culations in a model-independent manner, remarkably, 
there are configurations in the minimal supcrsymmetric 
extension of the standard model (MSSM) with the right 
properties for the effect we discuss here. In particular, we 
find that if the LSP is a neutralino quasi-degenerate in 
mass with the lightest stau, we can naturally obtain, at 
the same time, LDPs providing the correct dark matter 
abundance Q x h 2 = 0.113 [l| and NLDPs with the long 
lifetimes and the sizable densities in the early Universe 
needed in the proposed scenario. Such configurations 
arise even in minimal supersymmetric schemes, such as 
the minimal supergravity (mSUGRA) scenario plj and 
the minimal anomaly-mediated supersymmetry-breaking 
(mAMSB) model |22|. This implies that a detailed study 
of the (t,/^) parameter space using current and future 
cosmological data may constrain regions of the MSSM 
parameter space that are otherwise viable. Furthermore, 
we are able to make quantitative statements about test- 
ing the scenario we propose in future particle colliders or 
dark matter detection experiments. 

The paper is organized as follows: We first review in 
Section ^ how the standard calculation of linear per- 
turbations in an expanding universe must be modified 



to account for the effects of a decaying charged species, 
calculate the linear matter power spectrum and discuss 
the constraints to this model from big bang nucleosyn- 
thesis (BBN) and the spectrum of the cosmic microwave 
background (CMB). In Section lTTTl we briefly discuss how 
we estimate the nonlinear power spectrum from the lin- 
ear power spectrum and present several examples. In 
Section IIVI we discuss how measurements of the angu- 
lar power spectra of the primordial 21-cm radiation can 
be used to distinguish this effect from other modifica- 
tions to the primordial power spectrum. In Section 
we describe how this scenario can be embedded in a par- 
ticle physics model, concluding that the most appealing 
scheme is one where long lifetimes are obtained by con- 
sidering nearly degenerate LDP and NLDP masses. In 
Section IVTI we compute the lifetimes of charged next-to- 
lightest supersymmetric particles (NLSPs) decaying into 
neutralino LSPs, and show that, in the MSSM, the role 
of a NLDP with a long lifetime can be played by a stau 
only. In Section lVlIl we estimate what fraction of charged 
to neutral dark matter is expected in this case, while in 
Section IVIIII we describe how consistent realizations of 
this scenario can be found within the parameter space of 
mSUGRA and mAMSB. Finally, in Section|D[|we discuss 
the expected signatures of this scenario at future parti- 
cle colliders, such as the large hadron collider (LHC), 
and prospects for detection in experiments searching for 
WIMP dark matter. We conclude with a brief summary 
of our results in Section ^ 

II. CHARGED-PARTICLE DECAY 

In this section we discuss how the decay of a charged 
particle <j> (the NLDP) to a neutral particle \ (the LDP) 
results in the suppression of the linear matter power spec- 
trum on small scales. 

As the <fi particles decay to \ particles, their comoving 
energy density decays exponentially as 

p^a 3 = m<j > n ( p e~ t / r , (1) 

increasing the comoving energy density of x particles as 

p x a 3 = m x n X0 (l- Ue- t/T ). (2) 

Here n Xo — fl x p cr it/m x is the comoving number density 
of dark matter, and n^ Q = f$n Xo is the comoving number 
density of dark matter produced through the decay of cj) 
particles, a is the scale factor, and t is the cosmic time. 

Since the <f> particles are charged, they are tightly cou- 
pled to the ordinary baryons (the protons, helium nuclei, 
and electrons) through Coulomb scattering. It is there- 
fore possible describe the combined <f> and baryon fluids 
as a generalized baryon-like component (5 as far as per- 
turbation dynamics is concerned. We thus denote by 
Pp = Pb + Pip the total charged-particle energy density at 
any given time. At late times, after nearly all <f> particles 
have decayed, pp ~ pb- 
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FIG. 1: Shown is A 2 (fc) = k 3 P{k) /2tt 2 , the dimensionless 
matter power spectrum per logarithmic interval for the canon- 
ical n s = 1 ACDM model (dashed), for a model with r = 20 yr 
and f$ = 1/7 (long-dashed), for a model with r = 13 yr and 
ftf, = 1/5 (solid), and for a model with r = 3 yr and f$ = 1/2 
(dotted). Modes that enter the horizon be before the <fi par- 
ticles decay are suppressed by a factor of (1 — f^) 2 - 

The relevant species whose perturbation dynamics are 
modified from the standard case are the stable dark mat- 
ter (subscript x), the charged species (subscript /3), and 
the photons (subscript 7). By imposing covariant con- 
servation of the total stress-energy tensor, accounting for 
the Compton scattering between the electrons and the 
photons, and linearizing about a Friedmann-Robertson- 
Walker (FRW) Universe we arrive at the equations de- 
scribing the evolution of linear fluid perturbations of 
these components in an expanding Universe. In the syn- 
chronous gauge, for the '/?' component, the perturbation 
evolution equations are 

h = ( 3 ) 

and 

Op = --9p + c 2 s k 2 5p + ^Lan e a T (9 7 - 9 fj ) . (4) 

Here and in what follows, Sx — $Px/Px is the fractional 
overdensity and 9x — ikVx is the divergence of the bulk 
velocity in Fourier space of a given species X. An over- 
dot represents a derivative with respect to the conformal 
time. The number density of electrons is n e , while ot 
is the Thomson cross section. Because the <j) particles 



and the baryons share a common bulk velocity and over- 
density, 2 these equations are identical to the standard 
perturbation equations for the baryons with the replace- 
ment b — > (i (see, for example, Ref. [23^. For the dark 
matter we find that 

6 x = -9 x -h + \ m ^-(5 l3 -5 x ), (5) 
2 p x T 

and 

6 x = --e x + Xm E± ~(e -6 x ), (6) 

a p x r 

where A m = m x /m^ — (1 + Am/m x ) . The modifica- 
tions to the photon perturbation evolution are negligibly 
small because <j> decays during the radiation-dominated 
epoch when p^ -C p 1 and, as discussed below, for viable 
models X m ~ 1 to prevent unreasonably large spectral 
distortions to the CMB. 

Combining these equations with the (unmodified) 
equations for the neutrino perturbations we can solve for 
the linear power spectrum of matter fluctuations in this 
model. We have solved these equations using a modified 
version of cmbf ast j24[. In Fig. [I] we show the linear mat- 
ter power spectrum in this model for several values of the 
<f> lifetime r and fraction f$. As shown in this Figure, the 
small-scale density modes that enter the horizon before 
the <j> particles decay (when the age of the Universe is 
less than r) are suppressed relative to the standard case 
by a factor of (1 — /^) 2 - 

Since the decaying particles are charged, the produc- 
tion of the LDP will always be accompanied by an elec- 
tromagnetic cascade. The latter could in principle repro- 
cess the light elements produced during BBN, or induce 
unreasonably large spectral distortions to the CMB. We 
show here that in fact these effects are small for the mod- 
els discussed in this paper. 

The energy density released by the decay of <j> particles 
can be parametrized as 

(em — £EMf<t>Y x , (7) 

where e em is the average electromagnetic energy released 
in a decay and Y x = n x jn~ i is the dark-matter to 
photon ratio. In the specific models we discuss below, 
£em ~ Am/3, and 

Y x = = 3 x !„-" (™L) (She.) . (8) 

x m x n^ \m x J \0.23 J V ' 



2 This assumes adiabatic initial conditions. Note that the pertur- 
bation S^p = 5^ — 5p will generally evolve away from zero in 
an arbitrary gauge, even when starting with the adiabatic initial 
condition 5 C <j (0) = 0, due to gradients in the proper time. It 
is a special simplifying property of the synchronous gauge that 
Sab — for all time for adiabatic initial conditions. 
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This yields 

In the models we discuss below, Am/m x ~ 10~ 4 , and 
U < 1/2, giving Cba/ < 5 x 10~ 5 eV. The limit derived 
from too much reprocessing of the BBN light element 
abundances is (em < 3.8 Ty 1 / 4 x 10~ 3 eV [HIU where 
r yr = t/(1 yr), so we are safely below this bound. 

For r yr < 300, electromagnetic energy injection will 
result in a chemical-potential distortion to the CMB of 

m 

/Lt = 4.5 r y V 2 x 10- 3 e -°- 128 V 5/4 (10) 

and so we expect fi < 2.0 x 10 -7 — 2.3 x 10~ 6 for lifetimes 
between r vr = 1 — 100, below the current limit of /i < 
9 x 10" 5 111. 



III. THE NONLINEAR POWER SPECTRUM 

As density perturbations grow under the influence of 
gravity, linear evolution ceases to describe their growth 
and nonlinear effects must be taken into account. On 
large scales, where density perturbations have had insuf- 
ficient time to become nonlinear, the linear matter power 
spectrum describes the statistics of density fluctuations. 
However, on small scales the full nonlinear matter power 
spectrum is required. 

In order to calculate the nonlinear power spectrum 
for a given model we have used the recently devised 
halofit method |2^. This method uses higher-order 
perturbation theory in conjunction with the halo model 
of large-scale structure to determine the nonlinear power 
spectrum given a linear power spectrum. It has been 
shown to accurately reproduce the nonlinear power spec- 
tra of standard N-body simulations and, unlike the ear- 
lier mappings such as the Peacock and Dodds formula 
|30| , it is applicable in cases (like we consider here) when 
A 2 (k) = k 3 P(k) /2ir 2 is not a monotonic function. In par- 
ticular we have checked that it approximately reproduces 
the shapes of the nonlinear power spectra determined in 
Ref. [31j through N-body simulations in models where 
the linear power spectrum is completely cut off on small 
scales. As we are discussing here less drastic alterations 
to the linear power spectrum, we believe the halofit 
procedure provides an estimate of the nonlinear power 
spectrum adequate for illustrating the effect we describe 
in this paper. Any detailed study would require a full 
N-body simulation. 

In Figs. we show both the linear and nonlinear 
matter power spectra at redshift z = 4 (data from mea- 
surements the Lyman-a forest probe redshifts 2-4 at 
wavenumbers k/h ~ 0.1-10 Mpc -1 ) for a charged-decay 
model, and for a model with a running spectral index. 
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FIG. 2: Shown at a redshift z = 4 are A 2 (k) = k 3 P(k) /2tt 2 , 
the nonlinear (upper curves) and linear (lower curves) dimen- 
sionless matter power spectra per logarithmic interval for the 
canonical n 3 = 1 ACDM model (dashed), for an n 3 = 1 model 
with t = 20 yr and f$ = 1/7 (solid), and for a running-index 
model with n 3 — 1.00 and a 3 — —0.025 (dotted). Although 
the linear power spectra differ significantly in these latter two 
models, nonlinear evolution causes them to have nearly de- 
generate nonlinear power spectra for k/h > 1.5 Mpc - . 

Although these models have different linear power spec- 
tra, nonlinear gravitational evolution causes these models 
to have nearly identical nonlinear power spectra over in- 
teresting ranges of wavenumbers. The lifetimes shown 
were chosen because they produce effects on scales that 
can be probed by measurements of the Lyman-a forest, 
while the values = (1/2, 1/4, 1/5, 1/7) were chosen be- 
cause they arise in the supersymmetric models we discuss 
below. 

In Fig.[f)|we show the current constraints on the (n s ,a s ) 
parameter space from WMAP and SDSS Lyman-a for- 
est data. The constant-a s models we have compared 
our charged-decay model to lie within the WMAP+lya 
contours shown in the Figure. Since, in each case, the 
charged-decay model tends to interpolate between the 
standard ACDM model on larger scales (in particular 
on the scales probed by the CMB) and a constant-a s 
running-index model on smaller scales, and both of them 
are allowed by these data, we expect the charged-decay 
models we have considered in Figs. to be consistent 
with current data as well. We leave for future work the 
task of using a combined analysis of Lyman-a forest and 
other cosmological measurements to put limits directly 
on the (/^,t) parameter space (and thus on the parame- 
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FIG. 3: Shown at a redshift z = 4 are A 2 (k) = k 3 P{k) /2ix 2 , 
the nonlinear (upper curves) and linear (lower curves) dimen- 
sionless matter power spectra per logarithmic interval for the 
canonical n s = 1 ACDM model (dashed), for an n s = 1 model 
with t = 13 yr and = 1/5 (solid), and for a running-index 
model with n s = 1.00 and a a = —0.025 (dotted). Nonlinear 
evolution causes these two models to have overlapping non- 
linear power spectra for 2 Mpc -1 < k/h < 10 Mpc -1 . 



ter space of the MSSM models we discuss below). 



IV. THE 21 CM POWER SPECTRUM 

After recombination and the formation of neutral hy- 
drogen, the gas in the Universe cools with respect to the 
CMB temperature Tcmb starting at at a redshift z ~ 200. 
The spin temperature T s of the gas, which measures the 
relative populations of the hyperfine levels of the ground 
state of hydrogen separated by the 21-cm spin- flip transi- 
tion, remains collisionally coupled to the temperature 
of the baryons until a redshift z ~ 30 when collisions be- 
come inefficient and the spin temperature rises to Tomb- 
There is thus a window between z ~ 30-200 in which 
neutral hydrogen absorbs the CMB at a wavelength of 
21 cm. It has recently been suggested that the angular 
fluctuations in the brightness temperature of the 21-cm 
transition within this window may be measured with fu- 
ture observations and used to constrain the matter power 
spectrum at these very high redshifts pcj . 

At redshifts z > 30, the matter power spectrum on the 
scales k/h ~ 1 — 100 Mpc -1 of interest here are still in the 
linear regime. As discussed in Ref. j2jj, due to the un- 
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FIG. 4: Shown at a redshift z = 4 are A 2 (k) = k 3 P(k) /2tt 2 , 
the nonlinear (upper curves) and linear (lower curves) dimen- 
sionless matter power spectra per logarithmic interval for the 
canonical n 3 = 1 ACDM model (dashed), for an n 3 = 1 model 
with t = 15 yr and f$ = 1/4 (solid), and for a running-index 
model with n s = 0.96 and a a = —0.025 (dotted). Nonlinear 
evolution causes the charged-decay model to match the canon- 
ical ACDM model for k/h < 1.5 Mpc -1 and the running-index 
model for k/h > 1.5 Mpc -1 . 

precedented wealth of potential information contained in 
these 21-cm measurements, models with a running index 
or other small-scale modifications of the matter power 
spectrum can in principle be distinguished from each 
other. Even if the charged-particlc lifetime is so small 
that no significant modifications to the power spectrum 
occur on scales probed by other cosmological observa- 
tions (like the model shown in Fig. |3J, 21-cm observa- 
tions might detect or constrain such effects on the linear 
matter power spectrum. 



V. THE LONG LIVED CHARGED 
NEXT-TO-LIGHTEST DARK-MATTER 
PARTICLE 

From the particle-physics point of view, the setup we 
have introduced may seem ad hoc. We need a pair of 
particles that share a conserved quantum number and 
such that the lightest, the LDP, is neutral and stable, 
while the other, the NLDP, is coupled to the photon 
and quasi-stable, in order to significantly contribute to 
the cosmological energy density at an intermediate stage 
in the structure-formation process. Such a picture re- 
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FIG. 5: Shown at a redshift z = 4 are A 2 (fc) = fc 3 P(fc)/27r 2 , 
the nonlinear (upper curves) and linear (lower curves) dimen- 
sionless matter power spectra per logarithmic interval for the 
canonical n s = 1 ACDM model (dashed), for an n 3 = 1 model 
with r = 1 yr and = 1/2 (solid), and for a running-index 
model with n s — 0.96 and a s — —0.025 (dotted). Despite the 
drastic change in the linear power spectrum nonlinear evo- 
lution causes the charged-decay model to match the canoni- 
cal ACDM model for k/h < 8 Mpc" 1 and the running-index 
model for k/h>8 Mpc" 1 . 



quires three ingredients: (i) the relic abundance of the 
LDP must be compatible with the CDM component; (ii) 
the abundance of the NLDP must be at the correct level 
(namely, ~ 1/5 the total dark-matter density); and (Hi) 
the NLDP must have the proper lifetime (i.e., r ~ 10 yr). 

Let us start with the last requirement. One way to get 
the required lifetime is to introduce a framework with 
strongly-suppressed couplings. One such possibility is, 
for instance, to assume that the LDP is a stable super- 
weakly-interacting dark-matter particle, such as a grav- 
itino LSP in R-parity-conserving supersymmetric theo- 
ries or the Kaluza-Klein first excitation G 1 of the gravi- 
ton in the universal extra-dimension scenario |3^|. The 
NLDP can have non-zero electric charge, but at the same 
time a super-weak decay rate into the LDP (with the 
latter being the only allowed decay mode). In models 
of gauge-mediated supersymmetry breaking, this might 
indeed be the case with a stau NLSP decaying into a grav- 
itino LSP. In this specific example, we have checked that, 
to retrieve the very long lifetimes we introduced in our 
discussion, we would need to impose a small mass split- 
ting between the NLSP and LSP, as well as to raise the 
mass scale of the LSP up to about 100 TeV. This value 
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FIG. 6: We reproduce here (with permission) Fig. 3 from 
Ref. fl^| which illustrates the current constraints on the pa- 
rameter a s . The charged-decay models produce changes to 
the nonlinear power spectrum similar to constant a s ^ mod- 
els within the WMAP+Lya contour shown here and the ex- 
amples we show in Figs. |2HSl are denoted by the black circles 
(ours). However, the charged-decay models leave the CMB 
angular power spectra unaltered (corresponding to models 
with a s — on a line passing through the black cross). Match- 
ing to allowed constant-a s models is conservative in this re- 
spect, and we thus expect the charged-decay models we have 
discussed in Figs. |5J|S|to be consistent with current data. 



is most often considered uncomfortably large for a SUSY 
setup, and it also makes thermal production of LDPs or 
NLDPs unlikely, being near a scale at which the unitary 
bound [33] gets violated. Without thermal production 
one would need to invoke first a mechanism to wipe out 
the thermal components and then provide a viable non- 
thermal production scheme that fixes the right portion of 
LDPs versus NLDPs. Finally, such heavy and extremely 
weakly-interacting objects would evade any dark-matter 
detection experiment, and would certainly not be pro- 
duced at the forthcoming CERN Large Hadron Collider 
(LHC). This would then be a scheme that satisfies the 
three ingredients mentioned above, but that cannot be 
tested in any other way apart from cosmological obser- 
vations. 

An alternative (and to us, more appealing) scenario 
is one where long lifetimes are obtained by considering 
nearly degenerate LDP and NLDP masses. In this case, 
decay rates become small, without suppressed couplings, 
simply because the phase space allowed in the decay 
process gets sharply reduced. Sizeable couplings imply 
that, in the early Universe, the LDP and NLDP effi- 
ciently transform into each other through scattering from 
background particles. The small mass splitting, in turn, 
guarantees that the thermal-equilibrium number densi- 
ties of the two species are comparable. To describe the 
process of decoupling and find the thermal relic abun- 
dance of these species, the number densities of the two 
have to be traced simultaneously, with the NLDP, be- 
ing charged, playing the major role. This phenomenon 
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is usually dubbed as coannihilation |34j and has been 
studied at length, being ubiquitous in many frameworks 
embedding thermal-relic candidates for dark matter, in- 
cluding common SUSY schemes. 

We will show below that sufficiently long lifetimes may 
be indeed obtained in the minimal supersymmctric stan- 
dard model (MSSM), when the role of the NLDP is 
played by a stau nearly degenerate in mass with the 
lightest neutralino (with the former being the stable 
LSP and the thermal-relic CDM candidate we will fo- 
cus on in the remainder of our discussion). A setup of 
this kind appears naturally, e.g., in minimal supergravity 
(mSUGRA) [jjfj. This is the SUSY framework with the 
smallest possible parameter space, defined by only four 
continuous entries plus one sign, and hence also one of 
the most severely constrained by the requirement that 
the neutralino relic density matches the value from cos- 
mological observations. Neutralino-stau coannihilations 
determine one of the allowed regions, on the border with 
the region where the stau, which in the mSUGRA scheme 
is most often the lightest scalar SUSY particle, becomes 
lighter than the neutralino. Although a stau-neutralino 
mass degeneracy is not "generic" in such models, this 
scenario is economical in that the requirements of a long 
NLDP lifetime and of comparable LDP and NLDP relic 
abundances are both consequences of the mass degener- 
acy. In this sense, evidence for a running spectral in- 
dex or any of the other observational features we discuss 
would simply help us sort out which configuration (if any) 
Nature has chosen for SUSY dark matter. 



VI. LIFETIMES OF CHARGED NLSPS IN THE 
MSSM 

We refer to a MSSM setup in which the lightest neu- 
tralino Xi is the lightest SUSY particle. The charged 
particles that could play the role of the NLDP include: 
(1) scalar quarks, (2) scalar charged leptons, and (3) 
charginos. We now discriminate among these cases by 
the number of particles in the final states for the decay 
of NLSPs to neutralino LSPs. 

Scalar quarks and leptons have as their dominant de- 
cay mode a prompt two-body final state; i.e., S — > x^S 1 , 
where we have labeled S the SUSY scalar partner of 
the standard-model fermion S. A typical decay width 
for this process is 0(1) GeV, corresponding to a life- 
time 0(1CP 24 ) s. This holds whenever this final state 
is kinematically allowed; i.e., if mg > m$ + fn x o. If 
it is kinematically forbidden, there are two possibili- 
ties: squarks may either decay through CKM-suppressed 
flavor-changing processes or through four-body decays. 
For instance, the stop decay may proceed through t — > 
XiC or t — > Xibff- On the other hand, within the same 
minimal-flavor-violation framework, scalar leptons are 
not allowed to decay in flavor-changing two-body final 
states, and only the four-body decay option remains. 

The case for the chargino is different because this 



NLSP decay has a three-body final state, either with two 
quarks bound in a meson state — i.e., xt ~ * Xi 7r+ — or 
with a leptonic three body channel — i.e., xt - » xV- +v i- 
The latter final state becomes dominant, in particular, 
for electron-type leptons, I = e, if the mass splitting be- 
tween NLSP and LSP becomes small. 

We have listed all decay topologies as these are es- 
pecially relevant when discussing the limit in which we 
force a reduction of the allowed decay phase-space vol- 
ume; i.e., the limit in which the NLSP and LSP are quasi- 
degenerate in mass. Here we can also safely assume that 
the masses of the final-state particles, apart from the neu- 
tralino, are much smaller than the mass of the decaying 
particle. We can then consider the limit of a particle of 
mass m x o + Am decaying into a Xi ancl n — 1 massless 
final states, and derive an analytical approximation to 
the behavior of the final-state phase space d(f>^ and of 
the decay width r' n ) as functions of Am. In the case of 
two-body decays, the phase space reads 

d^> M fi-f ) 2 )ocAm. (11) 

Y 32?r 2 I ^m x o + AmJ v ' 

On the other hand, a recursive relation between d(f>^ 
and d^™ -1 - 1 based on the invariant mass of couples of 
final states yields 

/Am 
d^(d0 (2) )Gu) 

oc (Am) 2 ("- 2 ' +1 . (12) 

The dependence of the decay width T^ n > on Am must, 
however, take into account not only the phase-space de- 
pendence, but also the behavior of the amplitude squared 
MS n ' of the processes as a function of Am. The occur- 
rence of a massless final state yields, in the amplitude 
squared, a factor that scales linearly with the momenta 
circulating in the Feynman diagram. One therefore has 
the further factor 

M {n) oc (Am)"- 1 . (13) 

Finally, we have 

r (») K M (n) x d(j) (n) ^ (A. m ) 3 ™- 4 , (14) 

i.e., the lifetime to decay to a two-body final state scales 
like t^ 2 ) oc (Am) -2 , while for a four-body decay we have 
tW oc (Am)- 8 . 

Reducing the NLSP mass splitting Am one may hope 
to obtain "cosmologically relevant" NLSP lifetimes. For 
scalar quarks, this is not quite the case, as two-body final 
states always dominate, and even for amplitudes that are 
CKM-suppressed, the scaling with Am is too shallow, 
and the resulting lifetimes are rather short, even for very 
small mass splittings. 

The situation is slightly more favorable for charginos, 
whose three-body decay width is approximately equal to 
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to figure out a scenario in which the lightest smuon is 
lighter than the lightest stau. 

The scalings we have sketched are summarized in 
Fig. [7| where we plot lifetimes for a stop (CKM- 
suppressed), a chargino, and a stau as a function of Am. 
The decaying-particle masses have been set to 1 TeV, 
and the Am range is between the electron and the tau 
mass. The lifetimes of the stop and of the chargino have 
been computed with SDECAY [39j. Some details on the 
computation of the stau lifetime are given below. Notice 
that the scaling of Eq. (|14|l is accurately reproduced for 
all cases. The bottom line is that indeed the lightest stau 
can play the role of the NLDP with a cosmologically rel- 
evant lifetime, and that the stau is the only particle in 
the MSSM for which this can be guaranteed by adjusting 
only the LDP-NLDP mass splitting. 



FIG. 7: The lifetime of a 1 TeV stop (blue dotted line), a 
chargino (red dashed line) and a stau (black solid line), as a 
function of their mass splitting with the lightest SUSY parti- 
cle, the lightest neutralino. 
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At tree level and in the limit of pure Higgsino-like 
or Wino-like states, the lightest neutralino and light- 
est chargino are perfectly degenerate in mass. How- 
ever, when one takes into account loop corrections to 
the masses, Am usually turns out to be larger than a 
few tens of MeV; see, e.g., [33, |3^| . This translates into 
an absolute upper limit on the chargino lifetime of about 
1CH 2 - 3 ' s. 

Finally, in the case of sleptons I, the interesting regime 
is when Am < mi and only four-body decays are allowed. 
For example, for the lightest stau f{, the processes 



with 



fl - X>rff', 



(/',/) = ("n,fJ>), Oe,e), (u,d), (c,s), 



(17) 



where, depending on Am, only final states above the 
kinematic threshold are included. In this case, Am can 
be safely taken as a free parameter of the theory: in 
most scenarios the gaugino mass parameter setting the 
neutralino mass for bino-like neutralinos, and the scalar 
soft mass parameter setting the stau mass, are usually as- 
sumed to be independent; analogously, there are models 
in which the /j, parameter setting the mass for a Higgsino- 
like neutralino, and scalar soft mass parameter are un- 
related. A similar picture applies to the smuon, though 
lifetimes start to be enhanced at much smaller mass split- 
tings (m M instead of m T ), and it is theoretically difficult 



VII. THE RELATIVE ABUNDANCE OF THE 
CHARGED NLSP 



The next step is to determine the relic density of the 
NLSP and LSP in the early stages of the evolution of the 
Universe. As we have already mentioned, we are going 
to consider thermal production. We briefly review here 
how to compute the evolution of number densities with 
coannihilation |34L l35| . 

Consider a setup with N supersymmetric particles xi j 
X2j ■■■ Xn, each with mass mj and number of inter- 
nal degrees of freedom g,. The ordering is such that 
mi < m-2 < • • • < mjv- In the evolution equations, the 
processes that change the number density of SUSY par- 
ticles i are of three kinds: 



(a) Xi 

(b) X* 

(c) Xj 



> - X( , 

■ x l xi , 



Vj, 

V j ^ i, 

V j > i , 



(18) 



(16) 

where X a , X\, X(, and X[ are (sets of) standard-model 



L 6 ) 

(SM) particles. In practice, the relevant processes one 
should include are those for SM particles that are in ther- 
mal equilibrium. Assuming the distribution function for 
each particle k is the same as for the equilibrium distri- 
bution function, 



f k {E k ) cx fl\E k ) 



1 



exp(E k /T) ± 1 



(19) 



and invoking the principle of detailed balance, the Boltz- 
mann equation for the evolution of the number density 
of SUSY particle i, n, = <7i/(27r) 3 J d 3 pfi(E), normalized 
to the entropy density of the Universe, Yi = n,/s, as a 
function of the variable x = m\/T (with T the Universe 
temperature; this is equivalent to describing the evolu- 
tion in time) is given by 



In this equation, analogously to Y k , we have defined 
Y£ q = n cq / s 1 the ratio of the equilibrium number den- 
sity of species k (at temperature x) to the entropy 
density. On the left-hand side, we introduced g(x) = 
[1 + T/(3 g cS ) dgcs/dT}- 1 with g cS (T) being the effective 
degrees of freedom in the entropy density. The function 
g(x) is close to 1 except for temperatures at which a back- 
ground particle becomes nonrelativistic. On the right- 
hand side, the last two terms contain factors in T k ^i that 
label the partial decay width of a particle k in any final 
state containing the particle I; in our discussion they play 
a role just at late times when NLSPs decay into LSPs, 
giving the scaling we have used in Eq. |JT]|. The first two 
terms refer, respectively, to processes of the kind (a) and 
(b) in Eq. I|18|l , including all possible SM final and initial 
states. The symbol ((J a bV a b) indicates a thermal average 
of the cross section a a bVab] i-e., 

((JabVab) = ^ cq / d 3 p a d 3 p b f° q (E a ) f£ q (E b )(J ab V ab . 

As is evident from the form we wrote the Boltz- 
mann equation, interaction rates have to be compared 
with the expansion rate H of the Universe. In gen- 
eral, over a large range of intermediate temperatures, 
the (aix—ijVix—ij) Tlx terms will be larger than the 
{(TijUij) n eq terms, since we expect the cross sections to 
be of the same order in the two cases, but the scatter- 
ing rates will be more efficient as long as they involve 
light background particles X with relativistic equilibrium 
densities n c ^ that are much larger than the nonrelativis- 
tic Maxwell-Boltzmann-suppressed equilibrium densities 
n eq for the more massive particles j. This implies that 
collision processes go out of equilibrium at a smaller tem- 
perature, or later time, than pair-annihilation processes. 
Writing explicitly the expression for d/dx(Yi/Y^ q — 
Y k /Y^ q ), in which Maxwell-Boltzmann-suppressed terms 
and terms in mass splitting over mass scale can be ne- 
glected, one finds explicitly that in the limit that colli- 
sional rates are much larger than the expansion rate, for 
any i and k, Y i (x)/Y l cq (x) — Y k (x)/Y k eq (x), or equiva- 
lently, 

with Y(x) = Ek Y k(x) and Y cq (x) = £ fe ^° q (f)- At 
the temperature T c / , when Y^x( a *x->kVix^k) «x — H ■> 



collision processes decouple and the relative number den- 
sities become frozen to about 

n t (T) _ n q (T cfo ) _ g { ( m.A 3/2 / m k - mA 
n k (T) n cq (T cfo ) ~ g k \m k J \ T cfo J' 

(23) 

up to the time (temperature) at which heavier particles 
decay into lighter ones. 

The sum ^(a;) of the number densities has instead de- 
coupled long before. Eq. (|23|l is the relation that is imple- 
mented to find the usual Boltzmann equation [34| 35] for 
the sum over number densities of all species compared 
to the sum of equilibrium number densities, and that 
shows that the decoupling for Y occurs when the total 
effective annihilation rate becomes smaller than the ex- 
pansion rate, at a temperature T a /o that, as mentioned 
above, is much larger than T c f a . 

To get an estimate for T" c /o, we can take, whenever a 
channel is kinematically allowed, the (very) rough s-wave 
limit. 

. , 3 • lCT^cnr's- 1 

(Ti^kVi^k ~ cr lk v lk ~ {<r ik v lk ) ~ — — — , 

\l x h l 

(24) 

where an approximate relation between annihilation rate 
and relic density has been used [lfj. There are then 
two possibilities depending on whether (i) the back- 
ground particle X enforcing collisional equilibrium has 
a mass much larger than the mass splitting Am be- 
tween the SUSY particles involved, or (ii) the opposite 
regime holds. In the first case we find that T c f Q ~ 
mx/ (10 — 15), roughly the temperature at which X it- 
self (except for neutrinos) decouples from equilibrium. 
From Eq. I123H we find that riijn k ~ gi/g k ] i.e., they 
have comparable abundances. In the opposite case, 
we find instead that T c f ~ Am/(10 — 15) and hence 

/ s 3/2 

n-i/nk ^ gi/g k [y^j exp [(10 - 15)sign(m fc - m,-)]; 

i.e., the abundance of the heavier particle is totally neg- 
ligible compared to the lighter one. 

Long-lived stau NLSPs are kept in collisional equilib- 
rium with neutralinos by scattering on background 
and emission of a photon. In this case we are clearly in 
the limit (i), as m-pi — m x o <^ m T . Since the number of 
internal degrees of freedom for both staus and neutrali- 
nos is 2, we find that the fraction of charged dark matter 



in this model is 



9-n 



9* 



9r t 



(25) 



More carefully, though, this is not a strict prediction of 
our MSSM setup for a stau NLDP, and should be in- 
terpreted just as an upper limit. In fact, if we have 
other SUSY particles that are quasi-degenerate in mass 
with the LSP, and if they then coannihilate and decou- 
ple from the neutralino at a later time than the stau, 
either through mode (i) and then immediately decaying 
into neutralinos (in all explicit examples, the decay into 
staus is strongly suppressed compared to the decay into 
neutralinos), or through mode (ii), then is reduced to 



9r\ 



En 
i=i9i 



(26) 



where the sum in the denominator involves the neu- 
tralino, the stau and all SUSY particles with T c f a lower 
than the T c f for staus. 

In particular, if the lightest neutralino is a nearly-pure 
higgsino, then the next-to-lightest neutralino will also be 
a higgsino very nearly degenerate in mass, and the light- 
est chargino will also be nearly degenerate in mass, with 
mass splittings possibly smaller than m T . In this case, 
charginos and neutralinos will be kept in collisional equi- 
librium through scatterings on pairs. At the same 
time, the collisional decoupling of staus might be slightly 
delayed because of chargino-stau conversions through the 
emission of a photon and absorption of a tau neutrino; 
however this second process has a Yukawa suppression (as 
we are considering Higgsino-like charginos) compared to 
the first, and hence it is still guaranteed that the stau 
decoupling temperature is larger than the chargino T c f D 
temperature. Since the mass splitting between neutralino 
and chargino and that between lightest neutralino and 
next-to-lightest neutralino cannot be smaller than few 
tens of MeV (due to loop corrections to the masses), de- 
coupling will always happen in mode (ii) and we do not 
have to worry about possible stau production in their 
decays. Since g x + = 4, applying the formula in Eq. i|26|) 
we find = 1/5. In the case of a wino-like lightest 
neutralino, instead, the only extra coannihilating part- 
ner would be the lightest chargino, yielding = 1/4. 
Finally, adding to this picture, e.g., a quasi-degenerate 
smuon and selectron, = 1/7 could be obtained. 



VIII. LONG LIVED STAU NLSPS IN SAMPLE 
MINIMAL MODELS 

We provide a few examples of well motivated theoret- 
ical scenarios where neutralino-stau degeneracy occurs, 
possibly in connection with further coannihilating part- 
ners driving low values of f$. In surveying the possible 
models, the criterion we take here is that of the com- 
position of the lightest neutralino in terms of its dom- 
inating gauge eigenstate components. We thus outline 
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FIG. 8: Points of the mSUGRA parameter space, at null 
trilinear scalar coupling Ao, featuring f2 x o =0.11 and m x o = 
myj , for both negative (solid line) and positive (dashed line) 
sign of the fi parameter. 



bino-, higgsino- and wino-like lightest-neutralino bench- 
mark scenarios. 



A case with = 1/2: Binos in the mSUGRA 
model 



As we mentioned above, in the framework of minimal 
supergravity (mSUGRA) [21j one of the few cosmologi- 
cally allowed regions of parameter space is the tail where 
the neutralino and the stau are quasi-degenerate. In this 
case, coannihilations reduce the exceedingly large bino- 
like neutralino relic abundance to cosmologically accept- 
able values for neutralino masses up to around 600 GeV. 
Coannihilation effects depend on the relative mass split- 
ting between the two coannihilating species. Requiring a 
mass splitting as small as those found above amounts, as 
far as the neutralino relic density is concerned, to effec- 
tively setting rn x o = . This, in turns, sets the mass 
of the neutralino-stau system once a particular value of 
the relic abundance is required. We plot in Fig. [S] points 
fulfilling at once m x o = and f2 x o = JIcdm — 0.113, 
the latter being the central value as determined from the 



analysis of CMB data [l| . Result are shown in the (7 
tan 0) plane, with tan (3 the ratio of the vacuum expecta- 
tion values of the two neutral components of the SU (2) 
Higgs doublets, and at a fixed value of the trilinear cou- 
pling Aq = (this latter quantity is, however, not crucial 
here). The solid line corresponds to negative values of 
the Higgsino mass parameter /i, while the dashed line 
corresponds to positive values of fi. Notice that at fi < 
the accidental overlap of the heavy Higgs resonance with 
the coannihilation strip, around tan/3 = 43, shifts the 
neutralino masses to larger values. We point out that 
the two requirements of mass degeneracy and of the cor- 
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FIG. 9: A low-energy parameterization of the higgsino-like 
neutralino case, at tan/3 = 50. In the smaller frame we indi- 
cate the points, on the (u, Mi) plane, that give Qh 2 ~ 0.113. 
The corresponding models are reproduced in the larger frame, 
where we plot the relevant mass splittings of the chargino- 
neutralino system. 

rect relic density determine, at a given value of tan /3, the 
required mass of the neutralino-stau system, thus solving 
the residual mSUGRA parameter space degeneracy, and 
making the present framework testable and predictive. 

In the minimal configuration we have considered, only 
the lightest neutralino and the lightest stau are playing 
a role, hence = 1/2. However, since the mass splitting 
between the lightest stau and lightest smuon and selec- 
tron is rather small, assuming a slight departure from 
universality in the scalar sector, two additional quasi- 
degenerate scalar particles can be obtained and the frac- 
tion of charged dark matter reduced to — 1/4. 



B. A case with = 1/5: Higgsino-like neutralinos 

When the /i term is lighter than the gaugino masses 
Mi and Af2, the lightest neutralino gets dominated by 
the higgsino component. This situation occurs, again 
within the mSUGRA model, in the so-called hyperbolic 
branch/focus point (HB/FP) region 0,E1> where large 
values of the common soft breaking scalar mass too drive 
fi to low values. In this region, scalars are naturally 
heavy, at least in the minimal setup; however, the occur- 
rence of non-universalities in the scalar sector |42j | may 
significantly affect the sfermion mass pattern. In partic- 
ular, in a SUSY-GUT scenario, soft breaking sfermion 
masses get contributions from D-terms whenever the 
GUT gauge group is spontaneously broken with a reduc- 
tion of rank |43j| . Light staus may naturally occur, for 
instance when the weak hypercharge D-term dominates 
and features negative values. In this case the hierar- 
chy between diagonal entries in the soft supersymmetry- 



breaking scalar mass matrices is m 2 E <C rnfj D q l . The 
m 2 L term may also be lowered in presence of additional 
D-terms originating from the breaking of further U{\) 
symmetries. 

The relic neutralino abundance is here fixed by the 
interplay of multiple chargino-stau-neutralino coannihi- 
lations. In the limit of pure higgsino, the dynamics of 
these processes fixes the value of fi yielding a given relic 
neutralino abundance. On the other hand, a mixing with 
the bino component along the borders of the HB/FP 
region may entail a larger spread in the allowed mass 
range, affecting the \i higgsino fraction. We sketch the 
situation in Fig. |5J where we resort, for computational 
ease, to a low-energy parameterization of the above out- 
lined scenario. The smaller frame shows the points on 
the (/k, Mi) plane that produce the required amount of 
relic neutralinos. The larger frame reproduces the val- 
ues of the chargino-neutralino and neutralino-next-to- 
lightest-neutralino mass splitting; the lines end in the 
pure-higgsino regime. Suitable models, in the present 
framework, must also fulfill the mass splitting require- 
ment to- ~+ — to^q < to t . This enforces the allowed 

X2,Xl Al 

neutralino mass range, at tan/3 = 50, between 870 and 
880 GeV. Had we lowered the value of tan f3, the corre- 
sponding m x o range would only have shifted to masses 
just a few tens of GeV lighter. 

As we have already mentioned, since we are dealing 
with a case with two neutralinos, a chargino and a stau 
quasi-degenerate in mass, we find = 1/5. Again, a 
smuon and a selectron can be added to this to shift the 
charged particle fraction to = 1/7. 



C. A case with = 1/4: Wino-like neutralinos 

A benchmark case where the lightest neutralino is 
wino-like is instead provided by the minimal Anomaly 
mediated SUSY breaking (mAMSB) scenario j^. In this 
framework, tachyonic sfermion masses are cured by pos- 
tulating a common scalar mass term, too- The lightest 
sfermion turns out to correspond to the lightest stau f[ 
again. The latter, at suitably low mo values, may be 
degenerate in mass with the lightest neutralino. We per- 
formed a scan of the mAMSB parameter space, requir- 
ing to x o ~ TOfj , and found that the correct relic abun- 
dance requires the neutralino masses to lie in the range 
1250 < m x o < 1600 GeV, with the lower bound hold- 
ing for small values of tan (3, while the higher for larger 
ones. The chargino-neutralino mass splitting is always 
below the r mass, since \i m mAMSB is always a very 
pure wino, and it has large masses. We finally remark 
that here, as in the case of higgsinos, the occurrence 
of stau coannihilations raises the neutralino relic abun- 
dance, contrary to the standard result with a bino-likc 
LSP. 

Since, in this case, we have one neutralino, a chargino 
and a stau quasi-degenerate in mass, we find = 1/4. 






(a) 



FIG. 10: A few Feynman diagrams for four-body final states 
for the process t\ — * Xi^rff',- Diagram (a) is the dom- 
inant diagram; diagrams of type (b) are sub-dominant, and 
diagrams of type (c) are sub-sub-dominant. 



D. The Stau lifetime 

The t\ four-body decay proceeds through diagrams of 
the types sketched in Fig. QJJI They come in three sets: 
those with W exchange, those with exchange, and 
those with a sfermion exchange. However, since Am = 
m-fi — m x o is much smaller than any supersymmetric- 
particle mass, the virtuality of all diagrams except those 
featuring a r exchange (diagrams (a) and (b) in the 
Figure) is extremely large. Hence all diagrams but 
those with a r exchange will be suppressed by a fac- 
tor (m r /msuSY) 4 ~ 10~ 8 , and the interferences with the 
dominating diagrams by a factor (m T / Vtisusy) 2 ~ 10 ~ 4 . 
Of the two diagrams with a r exchange, however, the one 
with the exchange has a Yukawa suppressed H^ff' 
vertex, which gives a suppression, with respect to the 
VU-exchange diagram, 



(to m tan/3/mi4/) 2 (m T tan/3/m^/) 2 ~ 10 

for tan/3 



- 1CT 3 

5-50 
(27) 



in the most favorable muonic final channel. Notice that 
the chirality structure of the couplings entails that no 
interference between these two diagrams is present. Di- 
agrams with a charged Higgs or a sfermion exchange are 
moreover further suppressed with respect to those with a 
W exchange by a factor (mw / 'tosusy) which, depending 
on the SUSY spectrum, can also be relevant. 

In this respect, a very good approximation to the re- 
sulting stau lifetime is obtained by considering only the 
first diagram, whose squared amplitude reads 



E 



final states 



16(, 92 /V2) 4 (p„ T - Pf ) 
((p T ) 2 - m 2 T f {(p w ) 2 - mlyf 

\VR\ 2 ( 2 (Px? -Pr)(Pr -Pf>) ~ (Pr) 2 (p x » -Pp)j 

+m 2 .\V L \ 2 (p x a -pp) - 2m x om T Re(V L *V R )(p r -pp) 

(28) 

with Vl,r the left- and right-handed coupling of the T\ 
in the rirxi vertex. The sum is extended over the final 
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Wino, mAMSB Cm =1480 GeV) 
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10 100 
Mass Splitting (MeV) 



Electron 
Mass 



Tau 

Mass 



FIG. 11: The stau lifetime, as a function of the mass splitting 
with the lightest neutralino. The parameter space points are 



defined by the two requirements m 



and fL 



0.11. 



For the wino and higgsino cases tan /3 = 50, while in all cases 
Aq = and fi>0. 



states of Eq. (|17|) . For numerical purposes, the four- 
particles phase space is integrated with the use of the 
Monte Carlo routine Rambo, with final-state finite-mass 
corrections for all final states. 

Fig. ^2 shows the stau lifetime for a sample of the 
supersymmetric scenarios outlined in the preceding sec- 
tions. We fully account for threshold effects in the 
phase space, with the four contributions from electronic, 
muonic, and first- and second-generation quarks. The 
quark masses have been set to their central experimental 
values |44|. For the case of higgsinos, we reproduce the 
two extreme regimes when the m 2 L term is large (no Left- 
Right mixing) and when m 2 L ~ m? R (maximal Left-Right 
mixing). The differences in the lifetimes are traced back 
to overall mass effects and to the values of the Vl and 
Vr couplings (for instance, \Vl\ ^> \Vr\ in the wino case, 
while the opposite regime holds for the case of higgsinos 
and no LR-mixing). In any case, we conclude that life- 
times of the order of 1 — 100 years are obtained with a 
mass splitting Am = 20 - 70 MeV. 



IX. DARK MATTER SEARCHES AND 
COLLIDER SIGNATURES 



Unlike other charged long-lived NLDP scenarios, the 
framework we outlined above has the merit of being, in 
principle, detectable at dark-matter-detection and col- 
lider experiments. 

We show in Fig. 1121 the spin-independent neutralino- 
proton scattering cross section for the mSUGRA 
parameter-space points of Fig. [SI and for the higgsino 
and wino (mAMSB) cases, as discussed in sec. IVIII Bl 
and IVIII CI together with the current exclusion limits 
from the CDMS experiment ^^|, and the future reach 
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FIG. 12: The spin-independent neutralino-proton scattering 
cross section, as a function of the neutralino mass, along the 
parameter-space points outlined in Fig. |H| We also indicate 
the current exclusion limits from the CDMS experiment |45|. 
and the expected reach of the XENON 1-ton facility [i^l . 



of the XENON 1-ton facility gg. For binos, at fi > 
most points lie within less than one order of magnitude 
with respect to the future projected sensitivity, making 
it conceivable that this scenario may be tested in next- 
generation facilities. For negative /i, destructive interfer- 
ence among the lightest- and heavy-Higgs contributions 
lead instead to cancellations in a x o P . Finally, higgsino 
and wino detection rates respectively lie one and two or- 
ders of magnitude below the future expected sensitivity. 

Indirect-detection experiments look less promising, 
even for winos and higgsinos, which always feature un- 
comfortably large neutralino masses. We checked, for in- 
stance, that the expected muon flux from the Sun, gener- 
ated by neutralino pair annihilations, is at most 10~( 4 ~ 5 ) 
muons per km 2 per year, far below the sensitivity of fu- 
ture neutrino telescopes like IceCube g?} ■ 

Turning to collider experiments, considering as the 
searching tool the usual missing transverse-energy chan- 
nels, dedicated studies have shown that the mSUGRA 
coannihilation strip will be within LHC reach, mainly 
through in a mass range that extends up to m^o < 550 
GeV along the coannihilation strip, quite independently 
of tan fj [48| . Concerning higgsinos and winos, instead, 
the relevant mass range we study here appears to be be- 
yond standard LHC searches |57| . 

Even at high-energy colliders, the peculiar and distinc- 
tive feature of this scenario is however represented by 
the long-lived stau. The production of what have some- 
times been dubbed long-lived charged massive particles 
(CHAMP's) has been repeatedly addressed. Exclu- 
sion limits were determined by the CDF Collaboration 
[4^. l50| , and the future reach of the LHC and of a future 
Linear Collider has been also assessed for this broad class 



of exotic new particles [5lLl5^ |. The case of a stable 3 stau 
has been, in particular, considered several times, since it 
occurs in the context of various SUSY-breaking models, 
like gauge-mediated SUSY-breaking (GMSB) scenarios 
[53|]; see also, e.g., Ref. [3 EH- We remark that, con- 
trary to GMSB, in the present framework the production 
of long lived staus at accelerators is expected to come 
along together with the production of neutralinos, thus 
making the two scenarios distinguishable, at least in prin- 
ciple. 

A long-lived stau would behave as a highly penetrat- 
ing particle, appearing in the tracking and muon cham- 
bers of collider detectors with a small energy deposit in 
calorimeters. Staus, depending on their velocities, would 
produce either highly-ionizing tracks in the low-/3 regime, 
or, if quite relativistic, they would appear similar to en- 
ergetic muons. In this latter case, one could look at ex- 
cesses of dimuon or multi-lepton events as a result of su- 
perparticle production; for example, considering the ra- 
tio cr(/x + /i~)/(7(e + e~). In case additional particles have 
masses close to the neutralino-stau system, the total pro- 
duction cross section of superparticles would be greatly 
enhanced. Long-lived charginos may also give interesting 
accelerator signals |56| . 

While the reach of the Tevatron appears to be insuf- 
ficient to probe the parameter space of the models we 
are considering here [54J, the discovery of this kind of 
scenarios at the LHC, though challenging, looks quite 
conceivable, particularly in the case of binos. In fact, 
even in the less promising case in which the particle spec- 
trum does not feature any particle close in mass with the 
Xx — T\ system, the highly-ionizing-track channel should 
cover a mass range widely overlapping that indicated in 
Fig. [S] On the other hand, excess dimuon events could 
provide an independent confirmation, although the 5-cr 
LHC reach for CHAMP's in this channel alone has been 
assessed to lie around 300 GeV [54|. Moreover, if ad- 
ditional coannihilating particles (charginos, smuons, or 
selectrons) are present, the discover y a t the LHC would 
certainly look even more promising |55| . 

A further recently proposed detection technique for 
long-lived staus is represented by trapping these particles 
into large water tanks placed outside the LHC detectors 
[58| . Following the results of Ref. a 10-kton water 
tank may be capable of trapping more than 10 staus per 
year, if the mass ~ 400 GeV. The subsequent decays 
could then be studied in a background- free environment. 
More than twice as many sleptons would also get trapped 
in the LHC detectors, although in this case a study of the 
stau properties would look more challenging [58|. 



3 Here stable simply means that the decay length is much larger 
than the detector size. 
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X. CONCLUSIONS 

We have examined a scenario in which a fraction of 
the cold-dark-matter component in the Universe is gen- 
erated in the decay of long-lived charged particles and 
shown that a scale-dependent (or 'running') spectral in- 
dex is induced. The power spectrum on scales smaller 
than the horizon size when the age of the Universe is 
equal to the charged-particle lifetime gets suppressed by 
a factor (1 — /^) 2 . Such a feature might be singled out un- 
ambiguously by future measurements of the power spec- 
trum of neutral hydrogen through the 21-cm line, obtain- 
ing direct information on the charged particle lifetime 
and fy. 

On the contrary, current and future tests for depar- 
tures from a scale-invariant power spectrum based on 
Lyman-a data may fail to uniquely identify the scenario 
we propose. In fact, we have estimated the modifica- 
tions to the non-linear power spectrum at the redshifts 
and wavenumbers currently probed by Lyman-a forest 
data, and shown that these resemble (but are different 
in detail) those in models with a constant running of the 
spectral index a s . We expect, based on this resemblance, 
models with in the range 1/2 — 1/7 (as predicted in 
some explicit models we have constructed) to be compat- 
ible with current cosmological data for lifetimes in the 
range 1 — 20 yr. We have also verified that constraints 
from the primordial light-element abundances and dis- 
tortions to the CMB spectrum are not violated. 

From the particle-physics point of view, we have shown 
that the proposed scenario fits nicely in a picture in which 
the lightest neutralino in SUSY extensions of the stan- 
dard model appears as the cold-dark-matter candidate, 



and a stau nearly degenerate in mass with the neutralino 
as the long-lived charged counterpart. A small mass split- 
ting forces the stau to be quasi-stable, since the phase 
space allowed in the its decay process gets sharply re- 
duced. At the same time, it implies that neutralino and 
stau are strongly linked in the process of thermal decou- 
pling, with the charged species playing the major role. 
Owing to these coannihilation effects, the current neu- 
tralino thermal relic abundance is compatible with the 
value inferred from cosmological observatations and, at 
early times, the stau thermal relic component is at the 
correct level. 

We have described several explicit realizations of this 
idea in minimal supersymmetric frameworks, including 
the minimal supergravity scenario, namely the supersym- 
metric extension of the standard model with smallest pos- 
sible parameter space. We have pointed out that charged 
dark matter fraction from 1/2 to 1/7 (or even lower) can 
be obtained and that stau lifetimes larger than 1 yr are 
feasible. We have also shown that some of the models we 
have considered may be detected in future WIMP direct 
searches, and discussed the prospects of testing the most 
dramatic feature of the model we propose, i.e. the pro- 
duction of long-lived staus at future high energy particle 
colliders. 
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